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Abstract 

Using  a  three-dimensional  computational  model,  numerical  simulations  are  performed  to  investigate  the  performance  characteristics  of  proton 
exchange  membrane  fuel  cells  (PEMFCs)  incorporating  either  a  conventional  straight  gas  flow  channel  or  a  novel  wave-like  channel.  The  simulations 
focus  particularly  on  the  effect  of  the  wave-like  surface  on  the  gas  flow  characteristics,  the  temperature  distribution,  the  electrochemical  reaction 
efficiency  and  the  electrical  performance  of  the  PEMFCs  at  operating  temperatures  ranging  from  323  K  to  343  K.  The  numerical  results  reveal  that 
the  wave-like  surface  enhances  the  transport  of  the  reactant  gases  through  the  porous  layer,  improves  the  convective  heat  transfer  effect,  increases  the 
gas  flow  velocity,  and  yields  a  more  uniform  temperature  distribution.  As  a  result,  the  efficiency  of  the  catalytic  reaction  is  significantly  improved. 
Consequently,  compared  to  a  conventional  PEMFC,  the  PEMFC  with  a  wave-like  channel  yields  a  notably  higher  output  voltage  and  power  density. 
©  2007  Published  by  Elsevier  B.V. 
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1.  Introduction 

Since  their  introduction  by  NASA  in  the  1960s,  proton 
exchange  membrane  fuel  cells  (PEMFCs)  have  attracted  con¬ 
siderable  interest  within  industrial  and  academic  circles  as  a 
potential  power  source  for  a  range  of  mobile  applications. 
Characterized  by  a  compact  size,  an  efficient  performance,  a 
straightforward  design  and  operation,  and  a  high  degree  of  envi¬ 
ronmental  friendliness,  it  seems  likely  that  PEMFCs  will  one  day 
replace  gasoline  and  diesel  internal  combustion  engines  as  the 
method  of  choice  for  vehicular  applications.  The  literature  con¬ 
tains  many  experimental  and  numerical  investigations  into  the 
transport  phenomena  and  electrical  performance  of  PEMFCs. 
For  example,  Springer  et  al.  [1]  developed  a  one-dimensional 
isothermal  model  of  a  PEMFC  and  showed  that  the  water  diffu¬ 
sion  coefficient,  electro-osmotic  drag  coefficient,  water  sorption 
isothermals  and  membrane  conductivities  all  varied  as  a  function 
of  the  membrane  water  content. 
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In  a  PEMFC,  the  design  of  the  flow  delivery  system  has  a 
critical  effect  on  the  cell  performance.  Various  flow  channel 
designs  have  been  proposed,  including  straight,  serpentine  and 
inter-digitated.  Whilst  each  configuration  has  its  own  particu¬ 
lar  advantages,  the  inter-digitated  design  has  attracted  particular 
interest  due  to  its  efficiency  in  expelling  the  liquid  water  byprod¬ 
uct  from  the  porous  electrode.  For  example,  in  [2],  two  dead-end 
flow  channels  were  used  to  force  the  reactant  gases  through 
the  diffusion  layer  via  a  convection  effect  rather  than  diffusion. 
It  was  shown  that  the  resultant  shear  force  acting  on  the  fluid 
prompted  the  forcible  ejection  of  liquid  water  from  the  pores 
of  the  membrane.  Many  sophisticated  numerical  models  have 
been  developed  to  enable  the  influence  of  the  flow  field  config¬ 
uration  on  the  transport  phenomena  and  electrical  performance 
of  a  PEMFC  to  be  evaluated  [3-6]. 

In  conventional  PEMFCs  with  straight  gas  flow  channels, 
the  reactant  gases  diffuse  into  the  catalyst  layer  through  the 
gas  diffusion  layer.  However,  the  diffusion  process  is  inherently 
slow,  and  hence  the  electrical  performance  of  the  PEMFC  is 
inevitably  limited.  In  an  attempt  to  resolve  this  problem,  Kuo 
and  Chen  [7-9]  proposed  a  novel  PEMFC  design  in  which  the 
gas  flow  channel  had  a  wave-like  structure.  Using  a  3D  computa¬ 
tional  model,  the  authors  showed  that  this  unique  channel  design 
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Nomenclature 

C  mass  fraction 

Cp  quadratic  drag  factor 

D  mass  diffusion  (m2  s-1) 

F  Faraday’s  constant  (96,485  C  mol-1) 

ze  current  density  (A  m-2) 

im  ionic  current  density  (A  m-2) 

j  transfer  current  density  (A  m-3) 

kc  condensation  rate  constant  (Is-1) 

&e  evaporation  rate  constant  ( 1  (atm  s) - 1 ) 

&p  permeability  (m2) 

k^f  effective  ionic  conductivity  of  membrane 
(S  cm-1) 

elective  electric  conductivity  of  the  GDL 
(S  cm-1) 

M  molecular  weight  (kg  mol- 1 ) 

n  charge  number  of  the  sulfonic  acid  ions 

P  pressure  (atm) 

R  universal  gas  constant  (8.314  J  (mol  k)-1) 

S  source  term 

T  Temperature  (K) 

u  velocity  in  the  v-direction  (ms-1) 

v  velocity  in  the y-direction  (ms-1) 

V  operating  voltage  (V) 

w  velocity  in  the  z-direction  (ms-1) 

Zf  species  valence 

Greek  symbols 

a  transfer  coefficient  for  the  reaction 

ft  coefficient  of  thermal  expansion 

s  porosity 

X  water  content  of  membrane 

p  density  (kg  cm-2) 

o  ionic  conductivity  ( 1  ( £2  m) - 1 ) 

r  tortuosity 

v  viscosity  of  flow  (kg  (ms)-1) 

0  phase  potential  (V) 

?7act  activation  overpotential  (V) 

Subscripts 

a  anode 

c  cathode 

CL  catalyst  layer 

e  energy 

eff  effective 

GDL  gas  diffusion  layer 

h\  gas  flow  channel  height 

I12  gas  diffusion  channel  height 

in  inlet 

k  kt\ i  component  of  fuel  reactant 

m  membrane  phase 

mem  membrane 

out  outlet,  outer 


ref 

reference 

w 

wall 

X 

v-direction 

y 

y-direction 

z 

z-direction 

improved  the  uniformity  of  the  velocity  and  temperature  distri¬ 
butions  within  the  channel  and  reduced  the  intersection  angle 
between  the  dimensionless  velocity  vector  and  the  temperature 
gradient,  thereby  improving  the  heat  transfer  characteristics  of 
the  fuel  cell. 

The  current  study  performs  a  series  of  3D  numerical  simula¬ 
tions  to  investigate  the  influence  of  the  forced  convection  effect 
[10,1 1]  induced  in  the  wave-like  gas  flow  channel  on  the  veloc¬ 
ity,  temperature  and  gas  concentration  distributions  within  the 
PEMFC.  The  electrical  performance  of  a  PEMFC  with  a  wave¬ 
like  gas  flow  channel  is  then  compared  with  that  of  a  fuel  cell 
with  a  conventional  straight  channel  at  temperatures  of  323  K, 
333  K  and  343  K,  respectively.  The  validity  of  the  numerical 
model  is  confirmed  by  comparing  the  results  obtained  for  the 
electrical  performance  of  the  PEMFC  with  a  straight  channel 
with  the  experimental  results  presented  by  Wang  et  al.  [12]. 

2.  Simulation  model 

The  simulations  performed  in  this  study  are  based  on 
a  steady-state,  single-phase,  multi- species,  three-dimensional 
mass  transfer  model  of  the  PEMFC.  Fig.  1  presents  a  3D 
schematic  illustration  of  the  fuel  cell.  As  shown,  the  major  com¬ 
ponents  include:  (1)  two  wave-like  flow  channels,  (2)  two  carbon 
paper  gas  diffusion  layers  (GDLs),  (3)  two  catalyst  layers  and 
(4)  a  porous  membrane.  The  geometric  and  physical  parame¬ 
ters  applied  in  the  simulations  are  summarized  in  Table  1.  Note 
that  some  minor  parameters  (e.g.  the  internal  resistance,  the 
material  properties,  and  so  forth)  are  omitted  from  this  table, 
but  can  be  found  in  the  literatures  [13-15].  The  operating  pres¬ 
sure  and  temperature  are  assumed  to  be  1  atm  and  323  K,  333  K 
or  343  K,  respectively.  Regarding  the  reactant  gases,  the  anode 
is  supplied  with  humidified  hydrogen  with  a  mass  fraction  of 


Gas  channel  Bipolar  plate 


Fig.  1.  Schematic  representation  of  PEMFC. 
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Table  1 


Geometric  and  physical  parameters  employed  in  numerical  simulations 


Quantity 

Value 

Gas  channel  depth 

0.5  mm 

Gas  channel  width 

0.5  mm 

Gas  channel  length 

100  mm 

Gas  diffusion  layer  thickness 

0.3  mm 

Catalyst  thickness 

0.05  mm 

Membrane  thickness 

0.125  mm 

Porosity  of  gas  diffusion  layer 

0.4 

Porosity  of  catalyst  layer 

0.28 

Permeability  of  gas  diffusion  layer 

1.76  x  10“um2 

Permeability  of  catalyst  layer 

1.76  x  10-11  m2 

Permeability  of  membrane  layer 

1.18  x  10-18  m2 

Tortuosity  of  gas  diffusion  layer 

1.5 

Tortuosity  of  catalyst  layer 

1.5 

Electronic  conductivity  of  gas  diffusion  layer 

53  Sm-1 

Electronic  conductivity  of  catalyst  layer 

53  Sm-> 

Inlet  temperature 

323  K,  333  K  and  343  K 

Operation  pressure 

1  atm 

Anode  fuel 

H2  (76.5%)  and  H20  (23.5%) 

Cathode  fuel 

02  (21.5%),  N2  (77.2%)  and 
H20  (2.1%) 

Relative  humidity  of  the  anode 

100% 

70/30%  H2/H2O,  while  the  cathode  is  fed  with  saturated  oxy¬ 
gen  with  a  mass  fraction  of  21/79%  O2/N2.  The  N2  gas  serves 
as  a  diluent  and  is  assumed  to  be  inert. 

In  general,  the  straight  gas  flow  channels  used  in  a  conven¬ 
tional  PEMFC  are  designed  to  maximize  the  area  of  the  catalyst 
surface  exposed  to  the  reactant  gases  and  to  provide  a  route  for 
the  liquid  byproduct  of  the  catalytic  reaction  to  exit  the  fuel  cell. 
In  the  current  PEMFC  with  wave-like  anode  and  cathode  chan¬ 
nels,  the  reactant  gases  experience  a  mixed  convection  effect, 
i.e.  natural  convection  as  a  result  of  the  elevated  temperature  in 
the  gas  flow  channel  [16,17]  and  forced  convection  as  a  result 
of  the  wave-like  profile  of  the  channel  walls.  Consequently,  the 
wave-like  gas  flow  channels  not  only  maximize  the  reaction  area 
and  allow  the  water  byproduct  to  drain  from  the  fuel  cell,  but  also 
enhance  the  gas  velocity  in  the  vertical  direction  and  therefore 
increase  the  transfer  of  the  reactant  gases  to  the  catalyst  layer. 

In  performing  the  simulations,  the  following  assumptions  are 
made: 

1 .  The  axial  velocity  and  temperature  of  the  reactant  gases  are 
uniformly  distributed  at  the  inlets  to  both  gas  flow  channels. 

2.  The  reactant  gases  are  incompressible  and  enter  the  gas  flow 
channel  in  a  direction  normal  to  that  of  the  channel  cross- 
section.  The  gas  flow  in  both  channels  is  laminar  and  has  a 
Reynolds  number  of  less  than  200. 

3.  The  GDLs,  catalyst  layers  and  membrane  are  isotropic  and 
homogeneous,  and  have  a  high  permeability  and  a  uniform 
porosity. 

4.  The  electrochemical  reaction  is  governed  by  the 
Butler-Volmer  kinetic  equation. 

5.  The  water  byproduct  of  the  electrochemical  reaction  is  in  a 
vapor  state. 

6.  The  membrane  is  impervious  to  the  reactant  gases. 

7.  The  fuel  cell  geometry  is  periodic  in  the  v-axis  direction. 


Note  that  for  computational  efficiency,  the  model  is  divided 
into  seven  separate  layers,  namely  the  upper  wave-like  channel, 
the  anode  GDL,  the  anode  catalyst  layer  (CL),  the  membrane,  the 
cathode  CL,  the  cathode  GDL  and  the  lower  wave-like  channel. 
In  general,  the  heat  and  mass  transfer  characteristics  of  PEMFC 
gas  flow  channels  can  be  modeled  using  conventional  mass  con¬ 
servation,  Navier-Stokes,  and  energy  and  species  conservation 
equations.  The  basic  3D  gas  transport  equations  for  a  general 
PEMFC  can  be  expressed  as  follows: 


Continuity  equation: 


3  u  dv  3  w 

—  T  —  +  —  =  0. 
dx  3  y  3  z 

Momentum  equation: 


£eff 


3  u  3  u  3  u 

u - b  v - b  w  — 

dx  3 y  3  z 


3  2u 


3  2u 


3  2u 


SeffdP 

a? 


(1) 


(2) 


£eff 


dv  dv  dv 

u - b  v - b  w  — 

dx  3 y  dz 


d2v 


SeffdP  fd2v  d2v  „  . 

“  eff  W+d?  +  9?,1  +5t" 


(3) 


(  dw  dw  dw\ 

£eff  dP  ( d2w  3 2w  d2w\ 

=  “7^+U£eff  +  7z^)+s' 

•  Energy  equation: 


(  dT  3 T  3 T\ 

fetff  /X7  ipT  d2T\ 

“77  (a ?+97  +  9?J  +5e 

•  Species  conservation  equation: 


dCk  dCk  3  CA 


Afc.eff 


3  2Ck 


+ 


3  2Ck 


+ 


3  2Ck 


dx2  dy 2  dz 2 

•  Charge  conservation  equation: 


3  0e  #e  #e 

£eff  U - b  V - b  W - 

e  1  dx  dy  dz 


+  Sc. 


=  —  Sfj,. 


(4) 


(5) 


(6) 

(7) 


The  analytical  formulae  for  the  source  terms  Su ,  Sv,  Sw,  Se, 
Sc  and  in  Eqs.  (2)— (7)  are  derived  in  Table  2.  Note  that 
in  these  formulae,  the  parameters  £eff,  Cp,  kv  and  Zf  denote 
the  effective  porosity,  the  quadratic  drag  factor,  the  permeabil¬ 
ity  and  the  valence  of  the  species,  respectively.  Furthermore, 
Z\eff  =  DksTi  represents  the  effective  diffusion  coefficient  of 


Table  2 

Analytical  formulae  for  source  terms  in  governing  equations  (Eqs.  (2)-(7)) 
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the  Mi  component  of  the  reactant  fuel  [18-20].  In  the  PEMFC, 
the  generation/consumption  of  the  chemical  species  and  the 
charge  transfer  take  place  only  in  the  catalyst  layer.  Therefore, 
the  source  terms  in  Eqs.  (6)  and  (7)  can  be  implemented  in 
accordance  with  the  principles  of  electrochemical  kinetics,  i.e. 


*Sh2  =  - 
So2  =  - 


Sr2o  =  — 


Janode 

IF  ’ 

Jcathode 

4  F  ’ 

Jcathode 


2  F 


(8) 

(9) 

00) 


where  j  denotes  the  transfer  current  density  and  is  derived  from 
the  following  Butler- Volmer  kinetics  expressions: 


7a  —  7a, ref 


Ch2 


1/2 


CH2>ref 

exp  (— |^act) 


aa  F 

exp  (  —rj act 


(11) 


J  c  —  Ja,ref 


Co2 


a^F 

eXP  1  ~Rf^a 


^02,ref 

/  ac  \ 

-exp  [-  —  Friact ) 

where  rjact  is  the  surface  over  potential  and  is  defined  as 


(12) 


7?  act  —  <^a,c  tym  Eoc> 


(13) 


in  which  <pa  c  and  cpm  denote  the  potentials  of  the  carbon  and 
membrane  phases  in  the  catalyst  layer,  respectively,  and  Voc  is 
the  reference  open-circuit  potential  of  the  electrode. 

The  phase  potential  equation  for  the  potential  and  current 
profile  is  given  by  [21]: 


(14) 


where  0  is  the  phase  potential  function,  and  crm  is  the  membrane 
conductivity  and  is  expressed  as  [22] : 


crm(T)  =  cr™f  exp 


(15) 


where  cr^f  is  the  reference  conductivity  of  the  membrane  and  is 
given  by 

=  0.005 139*  -  0.00326,  (16) 


k= 


0.043+17.81a— 39.85a2  +  36.0a3 
14+  1.4(a-  1) 


for  0  <  a  <  1 
for  1  <  a  <  3 


(17) 


in  which  a  is  the  water  activity  and  is  defined  as 

th2o  P 

a  =  - . 

T^sat 


(IB) 


100 
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In  Eq.  (18),  the  saturation  pressure  varies  with  the  tempera¬ 
ture  and  can  be  determined  directly  from  thermodynamic  tables 
or  via  the  following  empirical  expression  [23]: 

p  _  2Q— 2.1794+0.029537— 9. 1837 x  10_5r2+l. 4454 xl0_7r3 


3.  Boundary  conditions 

The  governing  equations  for  the  current  PEMFC  model  have 
an  elliptic,  partial  differential  form.  Consequently,  boundary 
conditions  are  required  for  all  of  the  boundaries  in  the  compu¬ 
tational  domain.  Due  to  the  conjugated  nature  of  the  modeling 
problem,  the  gas  flow  channel  surfaces  are  included  within  the 
numerical  solution  domain  and  are  treated  as  a  fluid  with  partic¬ 
ular  properties. 

The  boundary  conditions  can  be  summarized  as  follows: 


1 .  Gas  flow  channel  inlet 

•  anode  inlet: 

U  =  Min,  T  =  7in,  V  =  0,  Ch2  =  (+,  in, 
C+O  =  <+,0.  m,  (20) 

•  cathode  inlet: 


T  —  ^in>  O  —  0,  Cq2  —  ^02,in’ 


Cn2  —  ^N2,in’ 

interface  between  gas  flow  and  GDL: 


3  u 

dy 

3  w 
dy 


3  u 

=£eff,GDL  — 
y=h\  dy 


y=h+ 


dv 

dy 


y=h  i 


dv 

=£  — 
dy 


=  £  - 


y=h  i 


dw 

dy 


y=h+ 


interface  between  GDL  and  catalyst  layer: 


£eff,GDL 


3  u 
dy 


=  £  - 


y=h  2 


du 

dy 


y=h\ 


dv 

£eff,GDL  — 
dy 

dw 

£eff,GDL  +— 

3  z 


dv 

=  £eff,CL  — 
y=h~2  ^ 


y=h\ 


y=h  2 


dw 

=  £eff,CL  — 

dz 


y=h+ 


(21) 

y=hf 

(22) 


(23) 


Lower  surfacecathode  gas  channel: 

o' 

II 

3 

II 

» 

II 

s 

(28) 

^surface  —  298  K, 

(29) 

Tin  >  7w. 

(30) 

4.  Numerical  method 

The  governing  equations  for  the  current  modeling  problem 
are  coupled  and  nonlinear.  Therefore,  the  equations  are  most 
conveniently  solved  using  some  form  of  numerical  method. 
In  the  current  study  the  continuous  cell  domain  is  discretized 
into  several  control  volumes  and  the  dependent  variables  are 
calculated  at  the  center  of  each  computational  cell  to  ensure 
that  the  conservation  properties  of  the  physical  quantities  are 
maintained.  The  solution  procedure  is  conducted  using  gen¬ 
eral  purpose  CFD  software  with  the  SIMPLEC  (semi-implicit 
method  for  pressure-linked  equations  consistent)  algorithm- 
based  solver  [24] .  The  solution  procedure  commences  by  solving 
the  velocity  fields  in  the  momentum  equation.  A  pressure  correc¬ 
tion  process  is  then  applied  to  take  account  of  mass  imbalance 
effects.  Finally,  the  concentration  equation  and  the  potential 
equations  are  solved.  The  simulations  assume  that  the  inlet 
section  of  the  channel  is  fully  developed  hydro-dynamically, 
and  hence  a  fully  developed  velocity  profile  for  rectangular 
ducts  is  imposed.  Additionally,  a  forced  convection  regime 
[25,26]  is  imposed  within  the  computational  domain  and  the 
Navier-Stokes  equations  are  solved  under  the  assumption  of 
laminar  flow.  Before  performing  the  simulations,  a  parametric 
study  was  conducted  to  identify  a  suitable  3D  grid  mesh  capa¬ 
ble  of  generating  accurate  representations  of  the  temperature 
and  velocity  gradients  near  the  walls  in  order  to  obtain  detailed 
insights  into  the  electrochemical  reaction  and  mass  transport 
phenomena  in  the  PEMFC.  The  results  suggested  the  use  of 
181  nodes  in  the  v-direction,  49  in  the  y-direction,  and  11  in 
the  z-direction.  A  grid  dependency  test  was  performed  to  con¬ 
firm  that  the  resulting  numerical  solution  was  independent  of  the 
grid  density.  In  the  subsequent  simulations,  the  iterative  com¬ 
putations  were  terminated  once  the  value  of  the  residues  fell  to 
less  than  10-6.  The  computations  were  performed  using  a  PC 
with  a  3.2  GHz  Intel  Pentium  4  CPU,  1  GB  DDR  RAM  and  the 
Windows  XP  operating  system. 

5.  Results  and  discussion 


•  interface  between  catalyst  layer  and  membrane: 


u  =  v  =  w  =  Ck  =  0. 

(24) 

Gas  flow  channel  outlet 

du  dv  dw  3  T 

dx  dx  dz  dx 

(25) 

Upper  surfaceanode  gas  channel: 

u  =  v  =  w  =  0, 

(26) 

^surface  =  298  K. 

(27) 

Fig.  2(a)  and  (b)  illustrates  the  distribution  of  the  vertical 
component  of  the  flow  velocity  in  the  straight  gas  flow  chan¬ 
nel  and  the  wave-like  gas  flow  channel,  respectively,  under 
operating  conditions  of:  (1)  a  cell  potential  of  0.6  V,  (2)  a  cath¬ 
ode  inlet  velocity  of  0.1  ms-1  and  (3)  an  anode  inlet  velocity 
of  0.05  ms-1.  Fig.  3  illustrates  the  distribution  of  the  axial 
flow  velocity  component  under  the  same  conditions.  Comparing 
Fig.  2(a)  and  (b),  it  is  clear  that  the  periodic  wave-like  structure 
of  the  gas  flow  channel  yields  a  significant  increase  in  the  verti¬ 
cal  velocity  component  of  the  gas  flow.  The  increased  supply  of 
reactant  gases  to  the  catalyst  layer  enhances  the  efficiency  of  the 
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(b).  Wave-like  geometry. 


0.16 

0.14 


Fig.  2.  Velocity  field  in  y-direction  of  gas  flow  channels. 


catalytic  reaction  and  therefore  improves  the  performance  of  the 
PEMFC.  Meanwhile,  comparing  Fig.  3(a)  and  (b),  it  can  be  seen 
that  the  wave-like  surface  introduces  a  nozzle-type  effect  which 
accelerates  the  flow  in  the  axial  direction  in  the  constricted  region 
of  the  channel  immediately  above  each  wave  feature.  However, 


U-m/s 

0.22716 


(a).  Straight  geometry. 

II  -  m/s 
0.41261 


(b).  Wave-like  geometry  . 


Fig.  3.  Velocity  field  in  x-direction  of  gas  flow  channels. 


GDL 


Flow  field  channel 


(a)  Wave-like  form  channel  T=323  K 

GDL 


(b)  Wave-like  form  channel  T=333  K 

GDL 


Flow  field  channel 

(c)  Wave-like  form  channel  T=343  K 


Fig.  4.  Liquid  water  saturation  held  in  cathode  gas  channel  and  diffusion  layer 
at  cell  voltage  of  0.6  V  and  inlet  temperatures  of  (a)  323  K,  (b)  333  K  and  (c) 
343  K. 

Fig.  3(b)  also  shows  that  the  axial  velocity  of  the  reactant  gases 
reduces  in  the  trough  regions  of  the  wavy  surface.  This  has  the 
effect  of  trapping  the  fuel,  thereby  increasing  its  supply  to  the 
catalyst  layer.  Figs.  2  and  3  both  show  that  a  strong  convective 
force  is  induced  along  the  reaction  surface  of  the  wave-like  gas 
flow  channel.  This  not  only  increases  the  supply  of  the  reactant 
gases  to  the  catalyst  layers,  but  also  improves  the  flow  of  the 
reaction  byproduct  out  of  the  PEMFC  and  therefore  reduces  the 
membrane  drowning  effect.  As  a  result,  the  performance  of  the 
fuel  cell  is  significantly  improved,  particularly  at  higher  current 
densities. 

The  liquid  water  saturation  field  in  the  gas  flow  channel  and 
the  GDF  on  the  cathode  side  of  the  PEMFC  at  inlet  cell  tempera¬ 
tures  of  323  K,  333  K  and  343  K,  respectively  are  as  Fig.  4(a)-(c). 
Comparing  the  three  figures,  it  is  observed  that  the  level  of  liquid 
water  saturation  increases  slightly  with  an  increasing  operating 
temperature.  Furthermore,  in  every  case,  it  can  be  seen  that  the 
liquid  water  saturation  increases  along  the  length  of  the  channel. 
As  a  result,  liquid  water  is  formed  at  the  interface  between  the 
membrane  and  the  cathode  catalyst  layer  in  the  outlet  region  of 
the  channel.  Consequently,  the  liquid  water  saturation  level  in 
the  catalyst  layer  is  higher  than  that  in  the  GDF.  Fiquid  water  is 
transported  from  the  catalyst  layer  towards  the  GDF  by  capillary 
action.  When  the  liquid  water  reaches  the  interface  between  the 
GDF  and  the  flow  channel,  it  is  transported  toward  the  channel 
outlet  by  the  drag  force  induced  by  the  convective  flow  of  the 
gas. 

Fig.  5  presents  the  variation  of  the  oxygen  concentration 
along  the  channel  in  both  the  wave-like  and  straight  gas  flow 
channels  at  a  cell  voltage  of  0.6  V  and  an  inlet  temperature  of 
323  K.  It  can  be  seen  that  the  oxygen  mole  fraction  decreases 
along  the  flow  direction  as  a  result  of  absorption  in  the  catalyst 
layer.  Furthermore,  the  results  confirm  that  a  greater  amount  of 
oxygen  is  consumed  in  the  wave-like  channel  due  to  the  forced 
convection  effect. 

Fig.  6(a)-(c)  illustrates  the  temperature  contours  within 
the  wave-like  gas  flow  channel  for  inlet  temperatures  of 
323  K,  333  K  and  343  K,  respectively.  The  results  reveal  that 
irrespective  of  the  inlet  temperature,  a  relatively  uniform  tem¬ 
perature  distribution  is  obtained  within  the  channel.  Hence, 
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Fig.  5.  Variation  of  oxygen  mole  fraction  in  cathode  catalyst  layer  at  cell  voltage 
of  0.6  V. 


Current  density  (A/cm2) 


Fig.  7.  Comparison  of  numerical  predictions  and  experimental  results  [12]  for 
variation  of  cell  voltage  with  current  density. 


it  can  be  inferred  that  the  wave-like  surface  is  instrumental 
in  improving  the  heat  transfer  within  the  channel.  In  gen¬ 
eral,  the  enhanced  convective  heat  transfer  performance  may 
be  the  result  of  the  increased  flow  interruption  caused  by  the 
wave-like  obstacles  in  the  channel,  a  reduction  in  the  thermal 
boundary  layer  in  the  region  of  the  channel  between  neigh¬ 
boring  waves,  or  an  increased  velocity  gradient  near  the  GDL 
boundary. 

In  this  study,  the  validity  of  the  numerical  model  was  verified 
by  comparing  the  predicted  results  for  the  polarization  charac¬ 
teristics  of  the  PEMFC  with  the  experimental  data  presented  by 
Wang  et  al.  [12] .  Fig.  7  illustrates  the  variation  of  the  cell  voltage 


(b)  Wave-like  form  channel  T=333  K 


Flow  field  channel 


(c)Wave-like  form  channel  T=343  K 


Flow  field  channel 


Fig.  6.  Temperature  contours  in  gas  flow  channel  at  cell  voltage  of  0.6  V  and 
inlet  temperatures  of  (a)  323  K,  (b)  333  K  and  (c)  343  K. 


with  the  current  density  for  both  the  wave-like  channel  and  the 
straight  channel.  Note  that  the  anode  and  cathode  channels  both 
have  an  inlet  temperature  of  323  K.  It  is  clear  that  a  good  agree¬ 
ment  exists  between  the  numerical  results  for  the  straight  gas 
flow  channel  and  the  experimental  data,  and  hence  the  validity 
of  the  computational  model  is  confirmed. 

Fig.  8  compares  the  polarization  curves  of  the  wave-like  chan¬ 
nel  and  the  straight  channel  at  inlet  temperatures  of  323  K,  333  K 
and  343  K,  respectively.  It  is  observed  that  a  good  agreement 
exists  between  the  two  sets  of  results  at  low  to  medium  current 
densities.  However,  at  high  current  densities,  the  wave-like  chan¬ 
nel  yields  an  effective  increase  in  both  the  cell  voltage  and  the 
power  density.  The  relatively  poorer  performance  of  the  straight 
gas  flow  channel  in  the  high  current  density  region  is  most  rea¬ 
sonably  attributed  to  the  lower  velocity  of  the  gas  fuel  supply  to 
the  GDF  and  catalyst  layers.  From  a  close  inspection,  it  can  be 
seen  that  both  the  cell  voltage  and  the  power  density  increase 
slightly  as  the  fuel  cell  temperature  increases.  An  improved 
PEMFC  performance  is  obtained  because  a  higher  temperature 
results  in  a  greater  catalytic  activity  and  a  greater  capacity  for 
water  removal  via  evaporation.  In  other  words,  increasing  the 
temperature  reduces  the  level  of  flooding  within  the  fuel  cell. 
Moreover,  a  higher  cell  temperature  increases  the  membrane 
conductivity  and  mass  diffusivity,  and  therefore  reduces  the  mass 
transport  resistance. 

Fig.  9  illustrates  the  variation  in  the  temperature  distribu¬ 
tion  in  the  axial  direction  of  the  straight  gas  flow  channel 
and  the  wave-like  gas  flow  channel  for  inlet  temperatures  of 
323  K,  333  K  and  343  K.  In  general,  it  is  observed  that  a  rea¬ 
sonable  agreement  is  observed  between  the  two  sets  of  results 
in  the  entrance  region  of  the  two  gas  flow  channels.  How¬ 
ever,  as  the  flow  distance  increases,  the  temperature  distribution 
within  the  wave-like  channel  remains  relatively  uniform,  while 
that  within  the  straight  channel  falls  rapidly.  The  relatively 
higher  temperature  within  the  wave-like  channel  enhances  the 
membrane  conductivity  and  accelerates  the  electrode  kinet¬ 
ics.  As  a  result,  the  electrical  performance  of  the  PEMFC  is 
improved. 
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Fig.  8.  Polarization  curves  as  function  of  fuel  cell  temperature. 


Fig.  9.  Variation  of  temperature  in  GDL  in  inlet  region  of  straight  gas  flow 
channel  and  wave-like  gas  flow  channel  for  cell  voltage  of  0.6  V. 

6.  Conclusions 

This  study  has  utilized  a  3D  computational  model  to  compare 
the  fluid  flow  phenomena  within  PEMFCs  containing  a  conven¬ 


tional  straight  gas  flow  channel  or  a  novel  wave-like  gas  flow 
channel.  The  simulations  have  considered  inlet  temperatures  of 
323  K,  333  K  and  343  K,  respectively,  and  have  focused  specif¬ 
ically  on  the  influence  of  the  forced  convection  effect  induced 
in  the  wave-like  channel  on  the  gas  flow  characteristics,  tem¬ 
perature  distribution,  electrochemical  reaction  efficiency  and 
electrical  performance  of  the  PEMFC.  In  general,  the  results 
have  shown  that  compared  to  a  straight  gas  flow  channel,  the 
wave-like  channel  provides  an  improved  convective  heat  transfer 
performance,  a  higher  gas  flow  velocity  and  a  more  uniform  tem¬ 
perature  distribution.  As  a  result,  the  efficiency  of  the  catalytic 
reaction  is  greatly  enhanced.  The  improved  catalytic  reaction 
results  in  a  higher  cell  voltage  and  an  improved  power  density, 
particularly  at  higher  values  of  the  inlet  temperature  and  current 
density. 
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